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Abstract. A quaternionic wavefunction consisting of real and scalar functions is 
found to satisfy the quaternionic momentum eigen value equation. Each of these 
components are found to satisfy a generalized wave equation of the form ^ d Q^±° — 
V 2 V>o + 2 (^) ^ + (^tt) 2 ipo = 0. This reduces to the massless Klein-Gordon 
equation, if we replace — > ^ + ■ For a plane wave solution the angular frequency 
is complex and is given by uj± = i m ° c ±cfc, where k is the propagation constant vector. 
This equation is in agreement with the Einstein energy-momentum formula. The spin 
of the particle is obtained from the interaction of the particle with the photon field. 
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1. Introduction 

Schrodinger equation has been generalized by Dirac and Klein-Gordon where in the 
former the resulting equation is a first order differential and in the latter is a second 
order differential equation in space and time. Dirac equation governs the motion of 
spin - 1/2 particles, while Klein-Gordon governs the motion of spin - particles. Dirac 
equation expresses the spin of the particle explicitly. Klein-Gordon equation doesn't 
explicitly reflect the spin of the particle. The spin of the particle shows up when a 
charged particle interacts with the electromagnetic field. In Schrodinger formalism 
the spin of the particle is obtained by replacing the momentum of the particle by its 
conserved conjugate momentum in the Hamiltonian of the system. Reduction of the 
Dirac equation for an electron in a magnetic field to its non-relativistic limit yields the 
Pauli equation with a correction term which represents the interaction of the electron's 
intrinsic magnetic moment with the magnetic field giving the correct energy. In this 
paper we would like to address the formulation of the Klein-Gordon equation. We 
have recently shown the richness of describing physical laws employing Quaternions, 
where we have found that the quaternionic Maxwell equations predict an existence of 
new scalar wave [1]. With the same fashion, we have unified the hydrodynamics and 
electromagnetism in an analogous way [2]. With the same enthusiasm, we would like 
here to write the quantum mechanical equations governing the evolution of microscopic 
world in terms of Quanternions. To this end, we use the eigen value equation, where we 
represent the wave function representing the microparticle, the energy and momentum 
operators as quaternions. It is interesting to note that the quaternionic wave function 
consists of scalar and vector components. The emerging equations can then describe 
scalar and vector particles. The resulting equations would then mix the scalar and 
vector components. We anticipate that these two particles could be the bosons and 
fermions. Hence, one can in some respect unify Klein-Gordon and Dirac equations in a 
single equation. The spin of the particle shows up when the particle interacts with the 
photon field. Moreover, the quaternionic uncertainty relation reduces to the ordinary 
uncertainty relations. This quaternionic relation yields the known Einstein's energy 
momentum relation. 

2. The model 

The quaternionic momentum eigenvalue problem can be written as 



(1) 



where 




(2) 



Equation (1) is in fact similar to Dirac equation [3] 

Pip = m cip , P = i ftrf^dp . 



(3) 



The Quaternionic Quantum Mechanics 



3 



Equation (1) is regarded as the quaternionic Dirac's equation. 

The product of the two quaternions, A = (a , a) and B — (b , bj is given by [4] 

a b — a -b , a b + ab + a x b) . (4) 



AB 

Apply Eq.(2) in (1) and use Eq.(4) to get 

V-^-^--^o = 0, (5) 

V^o-^-^ = 0, (6) 

and 

V x $ = , (7) 

where we have used the operators: p = —iKV and E = ih-^ as in ordinary quantum 
mechanics. Equations (5) - (7) are first order linear differential equations. They can 
be compatible with Lorentz transformation. The scalar equation above represents 
the energy equation, while the vector equation represents the momentum equation. 
In Dirac formalism, an electron, for instance, is described by a four component 
wavefunction (spinors). In the present formalism the particle is described by a 
quaternion wavefunction having also four components. Here, one component is a scalar 
and the three components are vector components. We may however ask the question 
that: are the two formalisms equivalent to each other? 

Equations (5) - (7) suggest that the wave functions; ip and ipo are not unique so that ip' 
and ip' can also satisfy Eqs.(5) - (6). This is possible if we define 
- OA 

^' = ^ + VA, ^=^ + _, (8) 

where A is some scalar function. Substituting Eq.(8) in Eq.(5) yields 

- . 1 d 2 A fm \ dk . . 

VA "sF-(f)aT = °- < 9 > 

if Eq.(5) is to be invariant under this transformation. Similarly, substituting Eq.(8) in 
Eq.(6) yields 

VA = 0, if m ^0, (10) 
so that Eqs.(9) and (10) imply that 

^+^A = 0, =► A = Aexp(-^t), A = const., (11) 
at n n 

hence 

^ = £, ^ = ^-^Aex V (-^ft). (12) 
However, if m = 0, then A satisfies the wave equation 

V 2 A = . 

c 2 dt 2 
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In this case the transformation in Eq. (8) is similar to the gauge transformation exhibited 
by the electromagnetic vector A and scalar <fi. Therefore, Eq.(5) - (7) are invariant under 
the transformation carried out in Eq.(12). 

Taking the V of both sides of Eq.(5) and using Eq.(6), (7) and the identity 
Vx(Vxf) = V(V • ■$) - V 2 ^, one gets 

^-^Mt)#+(T)^°- < 13 > 

Similarly taking the divergence of Eq.(6) and using Eq.(5), one obtains 

l^_W, + 2(^)^ + (^)\ = 0. (U) 

Alternatively, Eq.(14) can be derived if we differentiate Eq.(5) partially with respect 
to time and use Eq.(6). Equation (13) can be obtained if we differentiate Eq.(6) 
partially with respect to the time and use Eq.(5) and the vector identity V x (V x if;) = 
V(V ■ ip) — V 2 V>- Equations (5) and (6) represent vector and scalar waves. To my 
knowledge, they are new equations. It is interesting to see from Eqs.(14) and (15) that 
both wavefunction components satisfy the same wave equation. Moreover, Eqs.(13) and 
(14) are linear equations and hence they satisfy the superposition principle of quantum 
mechanics. For a massless particle (m = 0) they reduce to the ordinary wave equation, 
viz., V 2 ipo — Js-^ilr = 0. Notice that Eqs.(13) and (14) satisfy the transformation in 
Eq.(12). They can be written in operator form as 

irV=0, fl 2 V = 0, (15) 

where the operator fl is defined by 

c 2 dt 2 he dt h 2 y ' 

This can be casted in the form 

fl 2 = 4i-V 2 , where 9 (& ■ 



dr 2 dr \dt H 

It is thus apparent that Eqs.(13) and (14) when written in terms of the time r become 
Lorentz invariant. Equations (13) and (14) are similar to the evolution of the electric 
field in lossy medium [5]. 

If we consider a plane wave solution of the form 

ipo = C exp i(ut — k ■ r) , C = const., (17) 

and apply it in Eqs.(13) or (14), we obtain the dispersion relation 



2 o -m c 2 ( m c 4 2 2 
uj — 2% — - — cu — -— — h c k = . 



h \ h 2 

The solution of Eq.(18) is 



m ° c2 - Ml /-Ifrt 

u± = — - — i±ck. (19) 
n 
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This shows clearly that uj is complex. This is a quite interesting result. It can 
be interpreted as having two waves moving in opposite directions (because of the 
term, ±ck). The group velocity of these waves is the speed of light in vacuum 
(v g = = ±c). The phase velocity is given by v p = ^ = !Z ff-i ± c. The two 
solutions in (19) may also represent a particle with energy h~u + and antiparticle with 
energy hoj-. The difference of their energy states is AE = ±2h~kc = ±2h~u p , where u p 
is the photon frequency. This means that if one state is jumped to the other state two 
photons will be emitted (or absorbed). This interpretation resembles the Dirac theory 
of antiparticles. It is an interesting relation because it combines the wave and particle 
properties. It describes a matter wave associated with particles as in de Broglie theory. 
Now Eq.(19) yields 



ft>| 2 = m^c 4 + £rW, (20) 



so that 



m c = hk cS ,, k cS . = d — | k 2 , (21) 

is the particle's momentum in terms of angular frequency and wave number of the wave. 
Note that for light p = hk. Equation (20) represents the interaction between particle 
and wave properties. The refractive index of a medium is related to the frequency of 
the wave by the relation 

kc 1 m£ 
u = — n± = ± s--i ^—r, (22) 

where we used Eq.(19). The above equations shows that the refractive index is complex 
but |n| = 1. It is a rotation of the vacuum refractive index in the complex plane by an 



angle 



mpc 

a = tan" 1 y . (23) 

The imaginary term in Eq.(22) indicates the amount of absorption loss when the wave 
propagates through the space. 

Applying Eq.(19) in Eq.(13) and (14), one obtains 

Tfl C 2 

ipo = C exp( — t) exp i (± kct — k ■ r) , (24) 

h 

which represents a modulated wave with decaying amplitude. The physical significance 
of Eqs.(13) and (14) is still under investigation. The decaying factor in Eq.(24) depends 
on the momentum of the particle. It is apparent that fastly moving particle decays 
larger than slowly moving one. Hence, the motion of the particle distorts the wave 
associated with it. Therefore, a matter wave is a distorting light wave. Or equivalently, 
when light wave passes on a moving particle the wave is distorted by the motion of the 
particle. However, light wave is not distorted when encounters a stationary particle. 
Hence, Eqs.(13) and (14) represent the wave-particle interaction. It reflects the de 
Broglie hypothesis but in a rather different approach. 
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3. Klein-Gordon equation (KG) 

It is the equation of motion of a quantum scalar or pseudoscalar field, a field whose 
quanta are spinless particles. Any solution to the Dirac equation is automatically a 
solution to the Klein-Gordon equation, but the converse is not true. We would like in 
this section to deduce the KG from the set of equations given by (5) - (7). We seek the 
solution of the scalar part of the quaternionic wavefunction, viz., ip . We anticipate that 
the vector part, x/i, will govern the Dirac particles. Therefore, as in Dirac's equation 
case, in Eq.(14), we get 

This is a generalized Klein-Gordon equation with a time-dependent source term. The 
term on the right hand side of Eq.(25) can be seen as a dissipative (damping) term 
resulting from the inertia of the particle, or the motion of the particle in space-time. 
This implies that the inertia of the particle distorts the wave nature of the particle from 
a pure wave. For a massless particle Eq.(25) is a pure wave and the distortion term 
vanishes. It is evident from Eq.(25) that the damping interaction is inevitable for all 
non-zero mass particles. Now define the new time coordinate (r) as 

Applying this transformation in Eq.(25), we get 

a 2 V> = , where W 2 = -V 2 = U 2 . (27) 

c 2 or 2 

Thus, in the r coordinate the mass of the particle disappears. Or equivalently, due to 
the interaction (resistance) a massive particle appears to be massless. Thus, we can call 
the above transformation a killing mass transformation. Hence, it is interesting that 
a massive particle in one coordinate system can appear massless in another coordinate 
system. The electron is found to have an oscillatory motion with the extremely high 
frequency of uo = 2m ° c in addition to its normal linear motion. Schrodinger labeled 
that rapid oscillatory motion Zitterbewegung or "jittery motion" . Hence, the additional 
term in Eq.(26) may be attributed to this oscillatory nature of electrons. The oscillation 
is produced because of the interference between positive and negative energy states as 
described in Eq.(19). Several physicists proposed the idea that the electron spin and 
magnetic moment as generated by a localized circulatory motion of the electron. Using 
eq.(26), Eqs.(5) - (7) become 

V"0 - -5-5- = 0, 
c 2 or 

and 

V x $= 0. 
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Now write ipo(r,t) = ip{r) ip(t) and consider the wavefunction where ipo(r) = const.. 
In this case Eq.(25) reduces to a damped simple harmonic oscillator with a dissipative 
(damping) term proportional to the mass of the particle, i.e., 

dV + /m^V =Q (2g) 



or 



dt 2 h dt \ h 
This has a solution of the form 

2 2 

p(i) = Bexp(--^-t) + Dtexp(-^-t), : const., (29) 

<^(t) = Bexp(-yt) + J Dtexp(-yt), B , D : const., (30) 

The general solution of Eq.(25) as given by Eq.(24) now becomes 

ifj — Cexp(— ^ t) exp (± kct — k ■ r) i . (31) 

This wave depends on the rest mass of the particle and the propagation constant only. 
This shows that every mode with m 7^ damps exponentially to zero. 

We remark here that the energy equation in Eq.(20) coincides with the Einstein's 
relativistic energy-momentum formula, viz., E 2 = p 2 c 2 + m^c 4 , where hk = p. In this 
case Eq.(22) implies that a dissipative particle wave is described by the refractive index 

1 rape 

n± = ± — - % ^tt • (32) 



1 I '"0° 1 

1 1 k2l,2 1 - 



'()* 



h 2 k 2 x ^ h 2 k 2 

The refractive index approaches that of vacuum as the particle's mass approaches the 
speed of light in vacuum, i.e., n — > 1 as m — > 0. Alternatively, for the high energy 
particle ( h 2 k 2 » m c 2 ) one has n — > 1 as well. 

4. Interacting particles with electromagnetic radiation 

If the particle described by Eq.(2) interacts with an electromagnetic field defined by 

A= (%,/), (33) 



the momentum in Eq.(2) becomes 

P ->■ P + qA p -> p + qA and E ^ E + qip . (34) 

Applying Eq.(4) in Eq.(2) and using Eq.(33) yields 

V^- ? ^--^o = 0, (35) 

(36) 

Vx^=|-(W+i>o) , (37) 
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and 

tpip + c 2 A-if = 0. (38) 

These are two scalar equations and two vector equations. The scalar equation in Eq.(35) 
is the same as the scalar equation in Eq.(5). However, the vector equations are altered 
due to presence of the photon field. Now take the divergence of Eq.(35) to obtain 



c 2 dt 2 fx dt fx \ fx 2 fx 2 " c 2 

where we have used that fact that 



B = Vxf, V ■ (A x -0) = "0 ■ (V x A) - A ■ (V x ip) . 

The Klein-Gordon equation for a charged particle of mass rrx in an electromagnetic field, 
A^, reads 



2 2 \ 



D^ + -^-lVo = 0, where D„ = d„ + . (40) 

The covariant derivative in Eq.(40) can be compared with the transformation in Eq.(26) 
if we set qtp = m c 2 . Equation (39) can be written as 

1 d 2 ip 2 2m d^ qc T - / 'M c\ 2 

- V + -r — ?rr + — ^ ■ 5 + — - ^o = 0, (41) 



c 2 ^ 2 TU fx dt fx T \ fx 



where 



M 2 = m 2 + m 2 q , mq = ^A 2 -^, (42) 

is the effective mass of the interacting charged particle. The second term in Eq.(42) 
represents the mass of the photon felt by the charged particle due to the particle nature 
of the photon. It is interesting to see that this term is vanishingly small. For an electron 
with velocity v one has tp — v ■ A so that 



A v 2 

m q = q-^l-^. (43) 

For the photon (v = c), rrx q = 0. This is unlike the relativistic mass where rrx — > oo as 
v — > c. For a charged particle in a constant vector field (A), Eq.(41) reads 

1 d 2 ^ 2 2m o cdip , ( M c\ 2 

Upon using Eq.(26) the above equation reduces to 

rrx 2 c 2 

□ Vo + -^-^o = 0. 

This is a Klein-Gordon equation representing a charged massless (inertially) particle 
having an induced mass m q . Equation (42) represents the effective mass of the particle 
which is a function of space-time through A and </?. It is so amazing that the effective 
mass of such an interacting charged particle with electromagnetic field is directly 
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proportional to its charge. The fourth term on the LHS of Eq.(39) can also be linked 
to the spin of the particle, if we relate i/j to the spin S and i[)q by the relation 

f=^-S^o ^^-B= (±\ S-B^, (45) 



he h \h 2 . 

which corresponds to the interaction term (H int ) of the spin of the particle with the 
magnetic field of the photon, when a proper field theoretic model is found for the present 
formalism, viz., 

H in t. = (-)S - B = -il-B , fl=-^-S. (46) 

Using Eq.(45), the quaternion wave function can nw be defined as 

^=^o + ^=exp(-^S ? - J B)^ , (47) 

which means that the quaternionic wave function is a rotation (phase factor) of ipo- 
Applying Eq.(45) in Eq.(38) we get 

ip=~A-S. (48) 

Now upon using Eq.(26) in Eq.(41), we obtain 

a -> -> ui 2 c 2 

□'Vo + (^5 -B)^ + _^Vo = 0, (49) 

which represents the Klein-Gordon equation with a spin term. It is remarkable that the 
interacting Eqs.(35) - (38) reduce to the free Eqs.(5) - (7) when A = 0. 

Taking the divergence of E.q(36) and using Eqs.(35), (38) and the vector identity 
V • (fA) = /(V • A) + (V/) • A, one gets 

fev>-^)-f+fa-A-^)'k = 0. (50) 



dt ) T V c 2 dt 



If we replace A by -A, which is equivalent to choosing the complex field A* instead of 
A, the above equation yields 



E-tp = 0, (51) 



where 




E = ~[^ + ^\, V-l+i^ = 0. (52) 

Applying Eq.(51) in Eq.(45) implies that E ■ S = 0, i.e., the direction of the electric field 
is perpendicular to the spin of the particle. It is understood that the spin of the particle 
is connected with the magnetic field. However, there is no interaction of the particle 
spin with the electric field. Thus, Eq.(51) can account for this fact. Interestingly, Dirac 
just ignored this term, calling it unphysical for a point charge to have a polarization! It 
is very interesting here to see that Eq.(51) is satisfied if the Lorentz gauge is guaranteed. 
Hence, Eq.(49) is the generalized Klein-Gordon equation defining the spin of the particle. 
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This is unlike the ordinary Klein-Gordon equation where the spin of the particle is not 
accounted for. Applying Eq.(45) and (48) in Eq.(3) we obtain 

S x p = ihq ^A S 

This implies that S 2 = h 2 . This can result from the combination of two spin- 1/2 
particles. This supports the idea of supersymmetry that a boson results from two 
interaction fermions. Hence, the spin-orbit interaction term can be described by 

A-S 

(S x p ) = —L ■ S ' = ihq [ r ■ A ^— r ■ S 

where p, S and A are treated as operators. The average of the spin-orbit term 
vanishes. This implies that no such term can arise in a random motion of particles 
in an electromagnetic field. However, this can show up in an atom where it has been 
observed in Zeeman effect. 

5. The continuity equation 

The continuity equation governs the flow of any fluid. We would like to derive the 
continuity equation governing the particle flow represented Eqs.(5) - (7) and (35) - (38). 
Notice here that Eqs.(35) - (38) are the generalization of Eqs.(5) - (7). Multiply Eq.(35) 
from left by ipo and take the divergence of Eq.(36) with ip and add the resulting equations 
to get 

V-J + ^ = -^p, (53) 



where 



dt h 



f=ikj, p= ~(^ + \j\ 2 )=-^*%. (54) 



The RHS of Eq.(53) doesn't vanish as usual because of the dissipative effect of the 
mass of the particle. However, when m = 0, Eq.(53) reduces to the familiar continuity 
equation. It is interesting to note that the continuity equation is not affected by the 
presence of the fields (A, if). This merit is also exhibited by Brans-Dicke theory where 
a scalar field is introduced to couple to gravity in addition to graviton [6]. Now apply 
Eq.(45) in (54) to get 

, = _M( 1 + *V ( 55) 



he 2 c 2 \ h 2 t 

Since the current density is directly related to the spin of the particle, then this is a 
spin current density. In Eq.(55) p represents spin density. The current density and 
probability spin density are related by the familiar relation of if S = h, where we obtain 
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J = —pc. This implies that spin wave travels at speed of light. This equation can be 
written in the form 

-> -> dp 2m c 2 

V-J + ¥ = T = -^-p, (56) 

where T = — 2m ° c2 p is the torque density result from the two spin states. These two 
spin states behave like a dipole in the electric field of the photon. Hence, producing a 
polarization, P s given by T = —V ■ P s [7]. Hence, Eq.(56) is transformed into 

^ + V-(J + P s ) = 0, (57) 

which implies that the new current J e g. = J + P s is conserved. Alternatively, if we define 
a new time coordinate by the relation 
d d 2m c 2 

^ = ^ + — ' (58) 
then the continuity equation above reduces to 

^ + V.J = 0. (59) 
or) 

The factor of '2' appearing in the transformation in Eq.(58), when compared with 
Eq.(26), can be attributed to the fact that the particle acquires two energy states 
associated with its spin when it is placed in a photon field. This allow us to deduce that 
the particle has two spin-state. Hence, with such setting the continuity equation takes 
the normal form in the new time coordinate. 

6. Quaternionic uncertainty relation 

The quaternionic uncertainty relation can be written as 

(A?) (AP) > * . (60) 
Applying the multiplication rule in Eq.(4) yields 

h 

Af- Ap + AtAE > - , (61) 
2 

Af AE + c 2 At Ap = 0, (62) 

and 

AfxAp = 0. (63) 

Equation (62) implies that 

Ar 2 Ap 2 p 

At =C AE * V = C E> (64) 
where v — \-^\ is the velocity of the particle. Equation (64) is the combination of the 
two relativistic relations; E = mc 2 and p = mv, where m is the relativistic mass. It 
may also define the group velocity of the wavepacket, since 
d_E _ Ap_ _ J 

dp ~ AE ~ v [bb) 



The Quaternionic Quantum Mechanics 



12 



Equation (63) implies that 

Ay Ap z = Az Ap y , Az Ap x = Ax Ap z , Ax Ap y = Ay Ap x . (66) 

These equations imply that the fractional error in all components of the momenta and 
positions of the particle in any direction is the same (or the momentum gradient), i.e., 
l\x = ^Ky = ^ ne Q uari thy defines a current or mass rate. It thus implies 

that the mass current is conserved when a particle moves in space. It can be attributed 
to conservation of current or the continuity equation that can be written in the form 
V ■ p + ^ = for a particle with fixed volume. 

7. Conclusions 

We have investigated in this work the quaternionic eigen value problem in quantum 
mechanics. We have written the wave function, energy and momentum of the particle 
as quaternionic quantities. The eigen value equation reveals that the scalar and vector 
parts of the wavefunction are governed by a new wave equation. This is a new equation 
that we wish it will describe bosons and fermions. The quaternionic uncertainty 
relation predicts the Einstein relativistic formula of energy and momentum. When 
the interaction of the particle with an electromagnetic field is introduced a spin term 
appeared in the equation of motion. The scalar equation doesn't change due to the 
interaction of the electromagnetic field with the particle. However, the vector equations 
are altered. We generalized the ordinary uncertainty relation to a quaternionic one. 
This generalization provides us with the energy momentum relation of Einstein and 
the remaining uncertainty relations. Hence, by adopting the quaternionic quantum 
mechanics, namely, Dirac equation, we arrived at a dissipative or generalized Klein- 
Gordon equation with a particle spin. This formalism gives rise to the generation of 
spin angular momentum of the particle when a photon field is introduced in the equation 
of motion. Further investigation is going on to explore the physics of these two waves. 
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